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ON THE BREAKUP OF TECTONIC PLATES BY POLAR WANDERING 

by 

Han-Shou Liu 
ABSTRACT 

The observed boundary system of the major tectonic plates on the surface 
of the earth lends fresh support to the hypothesis of polar wandering. In this 
paper a dynamic model of the outer shell of the earth under the influence of polar 
shift is developed. The analysis falls into two parts: (1) deriving equations for 
stresses caused by polar shifting; (2) deducing the pattern according to which the 
fracture of the shell can be expected. For stress analysis, the theory of plates and 
shells is the dominant feature of this model. In order to determine the fracture 
pattern, the existence of a mathematical theorem of plasticity is recalled: it says 
that the plastic flow begins to occur when a function in terms of the differences of 
the three principal stresses surpasses a certain critical value. By introducing 
the figures for the geophysical constants, this model generates stresses which 
could produce an initial break in the lithosphere. The fracture pattern of plastic 
deformation in the outer shell of the earth due to shift of poles shows a remark- 
able correlation to the existing systems t f the tectonic plate boundaries. 
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ON THE BREAKUP OF TECTONIC PLACES BY POLAR WANDERING 
INTRODU CTION 

Convection currents in the interior of the earth (Pekoris, 1935) are capable 
of producing polar wandering, (Runcorn, Ib59). For the explanation of 
how a relatively thin shell has managed to slide a full 9«>° with respect to the 
underlying mantle, Goldreich und Toomre (1969) havt proposed a hypothesis 
that large angular displacement of the earth's poles L^s shifted on a geological 
time scale owing to the gradual redistribution of density inhomogeue'ties within 
the earth by the process of convection. Such a shift of polar axis with regard to 
the mantle must cause stresses and deformation of the thin shell because it would 
adjust its shape to the change of the flattening (Jeffreys, 1962). Therefore, the 
stresses and deformation generated by polar shift may provide a mechanism for 
the original breakup of the shell into several major tectonic plates. In the present 
paper, the equations for the stresses in a homogeneous shell of uniform thickness 
caused by a shift of the axis of rotation are derived. The magnitude of these 
stresses reaches a maximum value of the order of 10 dyn cm” which is suf- 
ficient for explaining a tectonic breakup. In order to deduce the fracture pattern 
according to which the brealcup of tectonic plates can be expected, we shall apply 
the theory of plastic deformation of shells (Hencky, 1924; Bijlaard, 1936; Van 
Iterson, 1943). The analysis of thi3 pattern gives an explanation of the existing 
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boundary systems of the major tectonic plates as described by Morgan (1968), 
LePichen (1968) and Isaoks, et al. (1968). 


STRESSES IN A SHELL 

Let r denote the radius of the earth, h the thickness of the shell, a and 
polar coordinates, a a , a Q and a r the normal components of stresses and r the 
shear stress. For an element prism as shown in Fig. 1, extending over the full 
thickness of the shell and bounded by two planes through the coordinate axis en- 
closing an angle da and two conic surfaces described around the axis with angle 
differing d/S, the conditions of equilibrium can be expressed by (Love, 1944; 
Timoshenko and Woinowsky-Krieger, 1959) 


<3cro 

a a hr cos fidfid a hr sin Bdfida. -o»h r cos fid/3 da 

o/3 ' p 

a) 

- ~ hr d/?da = 0, 
a a 


B<7 a 

da 


hr 


d/3da - 


— — hr sin /3d/3da - r hr cos /3d/3da 
“P 


- r hr cos fidfido. = 0, 


( 2 ) 


a f r 2 sin /3d/3da + hr sin /3d/3da + o a hr sin fidfida. - 0. (3) 
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The third terms appear in equations (1) and (2) because the lateral sides of 
the element shown in Fig. 1 have a trapezoidal form due to the curvature of the 
she'l. By simplifying these equations, we obtain 


( CT .. 


- C7g) COS 


Bcr„ 

fi - Sin ,3 + 


Br 

bI 


= 0, 


Ba 


Bt 


sin /3 - 2t cos fi - 0, 


( 4 ) 


cr , (cr a + = 0. 


The relations between stresses and displacements are 




T = 


2(1 


+ v)r\ 


Cot /3 


BU/j 


3«a \ 

+ ^' CSC/5 + U - 9 j ’ 




„ c sc/ 3-_ + U a c ot A 


3 )' 


( 5 ) 


in which E is Young's modulus, v is Poisson's ratio and U a , and U r are 
components of displacement in the elements of the shell. 

For the case of a change of the flattening, U a and r are zero in these co- 
ordinates and equations (4) and (5) reduce to 
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Bo - a 

— = (<^a " Ofl) Cot 

By 8 * * 


JS. 


h , v 

<? r = -— (^a + V* 


a ^' 


v ct„ 


E 

r \ B/3 



a a - V ~ (VjB Cot # + U r^* 


( 6 ) 


If we assume that the flattening of the earth's figure diminishes by an amount 
S and that the volume of the earth remains constant, U r can be formulated by 
(Jeffreys 1962) 


U r = Sr 


(• 


Cos 2 8 - 



(7) 


By substituting equation (7) into equation (6), the results are 

= ITTs ( 3 co,i 13 “ t)' 

u- = - 2(t/ j 12 8 - Sin /3 Cos A 
p v + 5 


( 8 ) 
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ft is noted that a t is small with respect to and o- 0 because the ratio 
4 h/r is about 6 x lO -2 . 

EFFECT OF POLAR SHIFT 

For deriving the stresses and deformation caused by a shift of the poles, we 
require the expression for the displacement U r in terms of latitude and polar 
motion before and after the event. As shown in Fig. 2, we consider the pole after 
the movement to point P at an angular distance 6 from its initial position P' . 

n n 

Denoting latitude by <j>, we may represent the radius of the earth by 

r = a (1 - A Sin 2 <P) • (®) 

If the earth has to adjust its shape entirely to the equilibrium figure after 
polar shift, the displacement U r becomes 


U r = Aa [Sin 2 <P - Sin 2 (< i> - 6)] 


( 10 ) 


= Aa Sin 6 Sin (2 <p - 9). 


U r reaches maximum value Aa sin 6 at <t > » n/4 + 6/2 and minimum value 
-Aa sin <9 at 0 = -rr/4 + e/2. The deformation of the outer shell due to polar 
shift can be obtained by the superposition of two special oases of deformations 
given by equation (7). The first case is the one with the axis coinciding with foe 
direction 0 » n/4 + 9/2 , and the second case is foe (me with foe axis coinciding 
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with the directUm <P - - -n /4 + 0 /2. Therefore, if we introduce successively 
/3 = <t> - (tt/ 4 + 0/2) and/3 * <t> - (-^/4 + 0 /2) in equation (7), we find 


8 


U r = Ar Sin 0 jcos* + §)] - Cos 2 [* - + f)]} 

= Ar Sin 0 Sin (20 - 0). 


( 11 ) 


which is identical with equation (10) to the first order of A. This result shows 
that we may derive the solution of our problem simply from equation (7). For 
this purpose we introduce a new system of polar coordinates 0 and K with its 
axis coinciding with the axis of rotation. Thus the origin of this system lies in 
the equator. We define </* from the great circle through the origin and the pole 
in its final position as shown in Fig. 3. According to equations (5) and (6), the 
position A t and A 2 in the original system of de f o rm ati o n have latitudes of 
w/4 + 0/2 and -tt/4 + 9/2 with regard to the initial position of the equator. 
Their latitudes with regard to its final position are tt/4 - 9/2 and -ir/4 - 9/2. 
Therefore, file coordinates of A t and A^ in the new system are 
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By introducing 8 = A sin/9 in equation (7) and (8), we obtain 

= \ frp + O = 7^-5 Sine ^ - Sin 2 ^ , 
= i (P B ~ O = AE Sine Sin 2 /3. 

£ V + J 


( 12 ) 


The general formula for stresses in a plane enclosing an angle y with that 
of the principal stresses are 


°X = °» + a d Cos 2 V » 


(13) 

^ = °- B - Otf Cos 2y, 
r = <7 a Sin 2 y. 


We denote /3, and y x for the system A x and /3 2 and y 2 for the system A 2 . 
The remits of o\,<^ and T for the combined system A t and A 2 in the new co- 
ordinates are 
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ca — AE Sin 0(Sin 2 /3. - Sin 2 /a.) 

* v + 5 

+ — L_ AE Sin 0(Sin 2 /3.Cos 2y.- Sin- (\ Cos 2y 2 ), 

v + 5 

y, = AE Sin 0(Sin 2 /0. - Sin 2 (i 2 ) (14) 

- — AE Sin 0(Sin 2 /3, Cos 2y. - Sin 2 & Cos 2y 2 ), 

v + 5 

r = — L_ AE Sin 0(Sin 2 & Sin 2y t - Sin 2 /3 2 Siny 2 ). 

2/ + S 

Since angle QOAj = (tt/ 4) - </* + (0/2) and angle QOA 2 = (3rr/4) - i/» + (0/2), 
we can derive 

Sin A Sin y, = Sin ^ - 'P + * 

Sin /3 a Sin y 2 = Sin - V' + = Sin + </< - 

(15) 

Sin A Cot yj = - Cos \ Cos - </» + » 

Sin A C*>* y 2 = - Cos \ Cos ^ - </» + ^ = Cos \ Cos ^- + 'Z' • 

By substituting equation (15) into equation (14), the results are 
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<7\ = — AE Sin 6 (Sin 2 k + 2) Sin (2i p - J) > 

^ v + 5 

CT = _1_ AE Sin 0 (3 Sin 2 *. _ 2) Sin (20 - 0), (16) 

w v + 5 

r = AE Sin # Cos K Cos (20 - #). 

v + 5 

Equations in (16) are >tresse*: in the outer shell of the earth caused Ly a 
shift of the poles over an angle 6 . 

BREAKUP PROBLEM 

Equation (16) shows that <?^ and vanish in two great circles through the 
origin given fay 0 - 6/2 and w/2 + @/2. Shear stress r in these circles 
reaches maximum value at the origin where the two circles intersect at right 
angles. The maximum value of t is 



For the shell of the earth* we adopt v - 0.25 and E » 10 12 dyn cm ~ 2 . 
(Gutenberg, 1959). At present the flattening of the earth A is 1/298. (Jeffreys* 
1962). Therefore* we have r mn * 1.8 x 10® dyn cm* 2 for 9 = (w/2). This shear 
stress Is quite considerable, and we may imagine that the outer shell of the earth 
will yield under this effect. 

For k » n/2* 0* n /4 + 6/2 and 0 * - 77/4 + d/2 * we conclude from equation 
(16) that r - 0. We fhrther find that - 3 &.<, everywhere in the meridian 
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K - (tt/ 2). The mairimnm absolute value o' ^ is found in the crossings of the 
meridian A. = (tt/2) with two great circles p = w/4 + 0/2 and p--tr/\ + 6/ 2. For 
the first case it is a maximum tension and for the second case a maximum com- 
pression. The absolute values in both cases amount to 


kl E .x = AESin5 - <*®> 

/v Btax 1/ + 5 

which gives |o\i aax = 1.9 x 10* dyn cm -2 for 6 = tt / 2 . This result also implies 
that the outer shell of the earth could not resist the stresses caused by a larger 
shift of the poles. 

At an early date of the earth's history the stresses generated by polar shift 
must have been even larger because the flattening must have been greater at that 
time. Using foe value A = 1/210 for a date 1.6 x 10® years ago (Jeffreys, 1962), 
we have * 1.8 x 10* <fyn cm" 2 and l°)J aM = 2.7 x lflP <fyn cm” 2 . There- 
fore, from foe estimation of foe stresses due to polar shifting, it is unlikely that 
the outer shell of the earth was aide to stand them without fracturing. 


ENERGY OF DEFORMATION 

The deformation energy which is required for the increase of the stresses . 
caused by polar shift can be expressed by the following formula (Timoshenko 
and Qoodter, 1961) 


V = 



J + “ jjg + o;) 2 • t 19 ) 


10 




13 


where <r t , cr y and a t axe components of normal stress, and r m , r y and r ( are 
components of shear stress in the rectangular coordinate system. For oar cue, 
equation (19) becomes 

<**> 

By sub s titut ing eqnatfon (16) into eqnatiop (20) and integrating oner the 
shell, we obtain 

V = . _L_ W A* Eh r 2 Sin* 6. (») 

IS v + 5 

DUfer en t ia t ing eq u at io n (21) wife resp e ct to #, the moment of farce retired 
to shiit the poles is 

T = - 7 ! • — irA a Eh r 2 Sin 2d. (2$ 

15 v + 5 

The maadnaun nhe of T occurs for e * w/4. Mrodnciag v = 046, 

E * 10 12 dyn cm* 2 , A = l/*98, h * 10 7 cm and r * x Iff cm, we fltad 
T b „ * 5.6 x 10 ,! 4ya cm. Therefore, the total drift moment of the convection 
o ar ren ts hi the upper mantle, as abutted by Iddye (1871) ad Bidder (1973), would 
be edeqaate to oseroome the elaatlc ooaaler foroa wtthlp tte ehell ac c ompa nying 
shift of the shell srowad fla earth. 
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GLOBAL FRACTURE PATTERN 

We consider now the problem of how the outer shell of the earth will deform 
and along which lines we may expect this to occur. For these deductions we apply 
the theory of plastic deformation of shells developed fay Hencky (1924), Bijlaard 
(1936), and Van fterson (1943). A c cording to this theory the plastic flow begins 
to occur when the stress ^ given by 

= ~ a 2> 2 + (°a - a $) 2 + < a i - • W 

surpasses a critical value ^ . In equation (23) <r t , <? 3 and a 3 are principal stresses. 
Bijlaard (1936), and Turcotte and Oxburgh (1973) have already applied the theory of 
plasticity for crustal problems and (harassed its validity in such cases. Intro- 
ducing foe stresses and r instead of ^ lt cr a and o 3 in equation (23), we find 

<r p - I + + 3t*| • (24) 

By iatrodacftng eqaatfton (16) in (24), the result is 

<T p = cr o tCOs a X + ^Sin (2^r - ff) Sin 4 \] 1/2 (25) 

where 

2' / 3AESin 6 
a s — * 

• V* 5 

Mace plastic deformation is characterised fay conditions of flow foe value 
of v in this case is about 0.5. 

U 
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Equation (25) describes the distribution of a p in the outer shell of the earth 
for each value of 6. It gives an idea about the zones of the shell that must first 
be liable to plastic flow. The value of it # is 2.1 x 10* dyn cm -2 for 6= 80*. ff 
we adopt toe critical stress er c = 1.6 x 10 9 dyn cm- 2 , the stress ratio I is 

1 = Jcos 2 \ + ^ Sin (2^ - 9 ) Sin 4 xj (26) 

Plastic fractur occurs when o- approaches a . Therefore, equation (26) 

P C 

describes toe plastic fracture lines in toe outer shell of toe earth caused by the 
shift of poles for I = 1. The fracture lines for 1=1 have been computed by means 
of equation (26) and have been drawn in stereographic projection as shown in 
Fig. 4. Fig. 4 can be changed over into another projection; we have transformed 
it onto a map in Fig. 5 for a shift of pedes over 80* along toe meridian of 75* W. 
This fr a cture pattern follows directions roughly in the N-E and N-W direction 
on toe northern hemisphere and in the S-E and S-W direction on tor <muthern 
hemisphere. Therefore the directions and locations of the fr ac t ur e lines in Fig. 

5 disclose the trend of the existing boundary system of six major tectonic plates 
(Eurasian, Indian, Pacific, North American, South American, African) as de- 
scribed by Morgan (19 68) LePichon (1968) and hacks et al. (1968). After toe 
original splitting Cf toe outer shell of toe earth as the first stage of sea floor 
spread ing (Blotter , 1973), the major tectonic plates poasihly became subjsct to further 
breakup by plate motions. Because the fracture lines are based upon the sim- 
plifying assumptions of toe uniform thickness and the homogeneity of the earth's 
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shell which are certainly not absolutely true, we have to expect that the boundary 
of tectonic {dates will not accurately follow the mathematical fracture lines which 
have been deduced. Taking this into account in the investigation of the cor- 
relation between the boundary features of plates and the fracture lines, we may 
consider them in good harmony. Since the choice of the magnitude and di- 
rection of the polar shift determines the position and direction of the f r actur e 
lines, the locations of plate boundaries at earlier stages in the earth's history 
may be investigated by following the path of polar wandering. 

CONCLUDING REMARKS 

A correlation has been found between the global boundary system of 
tectonic plates and polar wandering. The stress analysis in the present paper 
shows that the relationship of the fracture of the osier shell of the ear th with 
respect to the shift of pedes is influenced by geophysical parameters sach as the 
of the earth's figure* angnlsr shift of the axis of earth' s rotation* 
modulus of elasticity, mid Poisson's ratio of toe crust material of the earth, this 
and tee previous convection results, as ob t ained by Richter (1973) and IcMye (1971) 
can only have necessarily vague explanations of a tectonic c atastro phe owing to 
almost complete lack of under s t an ding of the breakup mechanism. They are* 
howevnr, not unexpected if plasticity arguments are accepted. Based on theories 
for tte of dofonutloii tadiitto wrrifOj, to rocuUiig fr&ctaro 

lines over ten surface of tee earth have been determined. These lines show a 
remarkable correlation to tfas boundary system of tee major tectonic {dates. If 
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this correlation is not fortuitous, and this does not appear probable, we have to 
suppose that the outer shell of the earth at some moment of its history has indeed 
u n dergone a corresponding plate-breaking in the process of polar wandering. 

The magnitude and direction of polar shift would alter die locations of frac- 
ture cones, ha this case the breaking mechanism of {dates would become dependent 
on the path of polar wandering. The problem of determination of the date for 
original breakup of the outer shell of the earth must therefore wait for detailed 
study of die evolution of piste motions. 

ACKNOWLEDGMENT 

1 thank R. Agrees for numerical analysis and computation. 

REFERENCES 

Btjlasrd, P. P n Theory of plastic deformation and its relation with negative 
gravity anomalies, oceanic trenches, geosyclines, volcaatem, orogency sad 
the West Pacific ocean. hoternai. Union of Geod. Geopfcys., Conference 
report, Edingbourg, 1998. 

Goidreich, P. and A. Toomre, Some remarks on polar wandering, J. Geopfays. 

Res., 74, 10, 2655-2567, 1969. 

G ut enbe r g, B n Physios of the Earth's interior, p. 180, Acad em i c Press, New 
York, 1959. 

Hsneky, H n Theory of plastic deformation, Proc. First Internet. Cong. App. 

Meolu, SIS, 1924. 

Mllyg . r/mHnn n ^nl Iwsthm hv MMteHnMMr Biafli OQCIfftOtiOB 

with d i fferen tial heating of tea orate, J. Qeophya. Bee., 75, 5, 1199-1158, 1971. 


15 


18 


Isacks* J. Oliver* and L. R. Sykes* Seismology and the new global tectonics* 

J. Geophys. Res., 73, 5855* 1968. 

Jeffreys, H n The Earth* Fourth edition* p. 328* Cambridge Univ. Press* I<ondon* 
1962. 

LePichon* X n Sea-floor spreading and continental drift, J. Geophys. Res., 73* 

3661. 1968. 

Love* A. E. H., A treatise on the mathematical theory of elasticity, p. 534* 

Dover Pub. Inc., Mew York, 1946. 

Morgan* W. J n Rise* trenches* great faults and crustal rocks, J. Geophys. Res.* 

73. 1959. 1968. 

Pekeris, C. L.* Thermal convection in the interior of the earth. Monthly Notices* 
Roy. Astron. Soc n Geophys. Supplement, 3* 8* 343-367* 1935. 

Richter* F. M.* Dynamical models for sea-floor spreading* Her. Geophys. Space 
Physics* 11* 2* 223-287* 1973. 

Runcorn* S. K n Rock magnetism* Science, 129* 1002-1012* 1959. 

Timoshenko* S n and S. Woinowsky-Krieger, Theory of plates and shells* p. 453, 
p. 533* McGraw-Hill Inc.* New York, 1959. 

Timoshenko, S. and .1. N. Goodier* Theory of elasticity, p. 148* McGraw-Hill 
Inc.* New York* 1951. 

Turcotte, D. L. and E. R. Oxburgh* Mid-plate Tectonics* Nature* 244* 337-339, 
1973. 

Van Iterson, F. K. Th„ Contribution to plasticity Theory* Ned. Akad* Wetenschap, 
52* 1* 1943. 


16 



. f f 1 1 

p - § § § 8 


19 


FIGURE CAPTIONS 

1. Stresses in spherical shell 

2. Effect of the shift of the earth's axis of rotation. 

3. The coordinate system \p and \ . 

4. Stress intensity I in stereographic projection of the half-sphere for 
ift of the poles over 9 . The solid lines for 1 = 1 are fracture lines. 

Figure 5. Fracture lines for a shift of the poles over 80° along the meridian 
of 75°W. The boundary system of six major tectonic plates coincides with 
the directions of these lines. 
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Figure 1. Stresses in spherical shell. 
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Figure 2. Effect of the shift of the earth’s axis of rotation. 
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system \p and 
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Figure 4* Stress intensity I in stereographic projection of the half-sphere for shift of the 
poles over 8 . The solid lines for I - 1 are fracture lines. 




Figure 5- Fracture lin^«. for u shift of the poles over 80 ° along the meridian of 7b W. The boundary system of six major 

tectonic plates coincides with the directions of these lines. 




